DE-FG05-97ER41031-57 



Hydrodynamics of the Vacuum 



P. M. Stevenson 

T. W. Bonner Laboratory, Department of Physics and Astronomy 
Rice University, P.O. Box 1892, Houston, TX 77251-1892, USA 



Abstract: 

Hydrodynamics is the appropriate "effective theory" for describing any fluid medium at 
sufficiently long length scales. This paper treats the vacuum as such a medium and derives 
the corresponding hydro dynamic equations. Unlike a normal medium the vacuum has no 
linear sound-wave regime; disturbances always "propagate" nonlinearly. For an "empty 
vacuum" the hydrodynamic equations are familiar ones (shallow water-wave equations) 
and they describe an experimentally observed phenomenon — the spreading of a clump 
of zero-temperature atoms into empty space. The "Higgs vacuum" case is much stranger; 
pressure and energy density, and hence time and space, exchange roles. The speed of sound 
is formally infinite, rather than zero as in the empty vacuum. Higher-derivative corrections 
to the vacuum hydrodynamic equations are also considered. In the empty-vacuum case the 
corrections are of quantum origin and the post-hydrodynamic description corresponds to 
the Gross-Pitaevskii equation. I conjecture the form of the post-hydrodynamic corrections 
in the Higgs case. In the 1+1-dimensional case the equations possess remarkable 'soliton' 
solutions and appear to constitute a new exactly integrable system. 



1 Introduction 



Quite generally, any fluid medium when viewed at sufficiently long length scales can be 
described by hydrodynamics. For length scales much, much greater than the mean free 
path the microscopic dynamics becomes irrelevant, except for determining the equation 
of state of the medium. The only relevant degrees of freedom are a density n and a flow 
velocity v: there is some 'stuff' and it flows. 

The idea of this paper is very simple: — to apply hydrodynamics to the vacuum. In fact, 
there are two cases; (i) an "empty" vacuum, where the equilibrium density is zero, and 
(ii) a "spontaneous condensate," or "Higgs-type" vacuum, where the relativistic energy 
density is minimized for some non-zero density n = n v . The actual physical vacuum 
is empty of certain particle species (e.g. electrons) but is a spontaneous condensate of 
others (e.g. quarks and gluons). The electroweak theory requires the vacuum to have a 
nonzero vacuum expectation value for the Higgs field, which one can view as a spontaneous 
condensate of scalar particles. This point has been well expressed by 't Hooft: "What we 
experience as empty space is nothing but the configuration of the Higgs field that has the 
lowest possible energy. If we move from field jargon to particle jargon, this means that 
empty space is actually filled with Higgs particles. They have Bose condensed." ^ |2] 

The great virtue of a hydrodynamic approach is that it is essentially independent of 
the microscopic dynamics. Hydrodynamics is almost as universal and fundamental as 
thermodynamics and it may yield important lessons for modern particle physics. (This 
view is advocated from a different perspective in an important recent paper In 
particular, hydrodynamics may perhaps provide a route to a deeper understanding of the 
"Higgs-type" vacuum, a phenomenon that is a vital, but experimentally untested feature 
of the Standard Model. 

Hydrodynamics, of course, has great limitations. It does not tell one how to produce 
or detect the excitations that it describes. Nor does it even tell one the scale of the 
phenomena, which is governed by a parameter whose value is set by the microscopic 
dynamics. I think it best to avoid speculations until the hydrodynamic equations and 
their solutions have been thoroughly explored. That task is begun in this paper, but 
much more remains to be done. For simplicity I concentrate on the 1+1-dimensional case, 
though a spherically symmetric d+ 1-dimensional solution will be discussed in Section 6. 
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2 Hydrodynamics of a normal medium 



In this subsection I very briefly review the basics of hydrodynamics applied to a "normal," 
nonrelativistic medium [IJ El |B] . I assume that the medium is "barotropic" (energy density 
only a function of the pressure). This is a good approximation for many normal media 
and will hold exactly in the vacuum case. As a further simplification I concentrate on the 
1 + 1-dimensional case. 

The hydrodynamic equations follow directly from conservation laws together with a 
constitutive relation for the medium, leading to coupled, first-order partial differential 
equations for p, the mass density, and v, the flow velocity. Mass conservation and momen- 
tum conservation yield the fundamental equations (subscripts indicate partial derivatives: 
pt = dp/dt, etc.): 

Pt + (pv) x = 0, (2.1) 

(pv) t + (pv 2 + P(p)) x = 0, (2.2) 

where P(p) is the pressure as a function of density. The previous equations together imply 
a simpler one, known as the Euler equation: 

v t + vv x + — = 0. (2.3) 
P 

Provided that the solution is everywhere smooth, one may regard Eqs. ()2.3|) and (|2.1|) 
as the basic equations. (However, if the solution develops discontinuities it is essential to 
remember that Eqs. (|2.1|) and (|2.2() are the fundamental pair.) 

For a normal medium the pressure varies linearly for small density disturbances: 

P(p) = P eq + v 2 (6p) + 0((6p) 2 ), (2.4) 

where 5p = p — p eq and Vq is a constant. For sufficiently small disturbances (5p <C p cq and 
v vq) the hydrodynamic equations can then be linearized, yielding 

Wt+«?^«0, (2.5) 
Peq 

(5p) t + PeqVx ~ 0. (2.6) 

If the x derivative of the first equation and the t derivative of the second are combined, 
so as to eliminate v, one obtains the wave equation for 5p: 

[5p) u - v%(5p) xx w 0. (2.7) 
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(Similarly, one can derive the wave equation for v.) Thus, in this linearized regime hy- 
drodynamics reduces to acoustics. Disturbances propagate as sound waves and obey the 
superposition principle. Note that the speed of sound is determined by the thermodynamic 
derivative 

(2.8) 

For larger disturbances, violating one or both of the conditions bp <C p e q, and v <C Vq, 
nonlinear effects come into play. In the vacuum case, as we shall see, there is no linear, 
acoustic regime; disturbances are always nonlinear. Insight into the nonlinear regime can 
be gained by considering first the case of a pressureless fluid, where the Euler equation 
reduces to Burgers' equation (also known by many other names) 

v t + vv x = 0. (2.9) 

The solution, easily verified by direct differentiation, can be expressed as follows [7j: Let 
the initial condition be v{x,t = 0) = V(x), then at later times 

v = V( X ) with X = x -V(x)t, (2.10) 

where the latter equation defines \ implicitly. Thus, velocity disturbances propagate 
locally at the local flow velocity. If V is monotonic increasing to the right then the forward 
part of the wave outruns the backward part and the disturbance simply stretches out. 
However, if V is monotonic decreasing then the backward part of the wave tends to catch 
up with the forward part, steepening the wave profile, as with ocean waves approaching 
the shore. After a finite time the solution, rather than becoming multivalued, as with 
waves at a beach, develops a finite discontinuity in v, a shock wave. 

Quite generally, the solution may develop shocks (finite discontinuities in p and/or v) 
in finite time, even from very smooth initial conditions. Requiring mass and momentum 
conservation yields two equations: 

, , , [Lewi] \[pv 2 + P(p)}} , s 

shock speed = ^rr = ,7 „ , ( 2 - n ) 

[[p}\ [Hi 

where [[X]] = A^ght — Ai e ft. The shock solution, even though the first derivatives do not 
exist at the shock, can be regarded as a "weak solution," a solution in a distribution- 
theoretic sense, of the original equations |S]. An important warning here is that the 
mathematical notion of a "weak solution" must be applied to the fundamental mass and 
momentum conservation equations, (|2.1|) and (|2.2|) . and not to the Euler equation (|2.3|) . 
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At a shock the solution violates the basic hydrodynamic assumption that all relevant 
length scales are much greater than the mean free path. Thus, corrections to hydrody- 
namics play a role. These corrections can be developed as a series expansion in (mean free 
path) / (length scale) , known as the Chapman-Enskog expansion [S] . For a normal medium 
the post-hydrodynamic corrections involve dissipative effects such as viscosity. The correc- 
tion terms involve second derivatives and violate (macroscopic) time-reversal invariance, 
since they entail entropy generation. In the pressureless-fluid case one encounters the 
viscous Burgers' equation [TU] . 



in which the viscosity term ev xx smooths out the shock discontinuity. For this equation, 
and more generally, one may find the shock profile by looking for "travelling wave" solu- 
tions with v = v(x — at) and p = p(x — at) where a = is a constant and vl ^vr and/or 
PL 7^ PR-, where vl =v(x-*- oo), etc.. As e — ► the travelling wave's profile narrows, ap- 
proaching a discontinuity, and the speed a tends to the shock speed of Eq. (|2.11|) . In this 
sense, a hydrodynamic solution with shocks is still a meaningful approximate description 
of the physics. 

3 Constitutive relation of the vacuum medium 

Henceforth a relativistic framework will be adopted; energy and energy density will include 
the rest-energy (mc 2 ) contribution and the speed of light, c, will be set to 1. 

The "empty vacuum" and "Higgs vacuum" cases are illustrated in Figs. 1 and 2. In 
both cases, since we are dealing with a zero-temperature medium, the pressure is given by 



where £ is the energy density and n is the number density of particles. Note that a precise 
definition of n is not needed. What matters physically is the "constitutive relation" - 
how P varies with £ in the vicinity of the equilibrium state, and it is merely a convenience 
to express this relation parametrically as £ = £(n) and P = P(n) via the above equation. 
If, for example, the underlying microscopic theory is a scalar quantum field theory, one 
may take n to be proportional to (j) 2 , where cj) is the classical field, and identify £ (n) with 
the field-theoretic effective potential V e s((p). 

The empty-vacuum case is essentially nonrelativistic. The energy density £ is domi- 
nated by the mass density (times c 2 , but c = 1 here). However, this linear term exactly 



v t +vv x = ev xx , 



(2.12) 



P = - 




(3.13) 
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Figure 1: The empty vacuum case; energy density as a function of mass density p = ran. 
The equilibrium density is p = 0. 
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C is a constant. 
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cancels in the pressure (j3.13|) . so that the small non- linear term \gp 2 is crucial jllj : 



£ = p + \gp 2 + ..., 



P 



The speed of sound is formally given by 



\9P 2 + 



dP 
~d£ 



p=0 



dP 
dp 



0. 



(3.14) 
(3.15) 

(3.16) 



p=0 



In the Higgs- vacuum case we find a "mirror-image" situation with the forms of £ and 
P interchanged. (For convenience I subtract a constant from the energy density so that £ 
vanishes in equilibrium. Since gravity is not being included in the discussion this step is 
quite innocuous.) 

£ = \Ca 2 + . . . , 



with [T2] 



a 



a + \Ca l + 



C 



71, 



where C is a positive constant. C is the compressibility of the vacuum medium 

\ -l 



C 



n 



,f£_ 

"dn 2 



(3.17) 
(3.18) 

(3.19) 
(3.20) 



For small disturbances one has P>£, so the Higgs- vacuum case is ultrarelativistic. The 
speed of sound is formally infinite (14} : 



dP 

~di 



dP/da 



(T=0 



oo. 



(3.21) 



a=0 d£ / d(J 

This is the first sign that something very interesting — not to say bizarre — is going on. 
Note that the vanishing denominator arises because d£/dn vanishes at n v , which is the 
defining property of a Higgs-type vacuum. 

To anticipate (and oversimplify), the empty- vacuum case behaves like Burgers equation 
Vt + vv x = and effects tend to "propagate" at the local flow velocity v , while the Higgs- 
vacuum case behaves like the x «-> t Burgers equation vt + -v x = whose implicit solution 
v = v(x — -t) implies that disturbances tend to "propagate" at speed A (that is c 2 /v, 
restoring c). Thus, the "propagation" speed is superluminal, and can be arbitrarily large 
as v becomes arbitrarily small. This is a fundamental and inescapable characteristic of 
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the Higgs-type vacuum. Obviously, deep issues related to causality are involved. Note 
that it is well known that an ultrarelativistic medium (where the pressure is much greater 
than the energy density) can have a superluminal speed of sound Discussions of the 
causality issue for an ultrarelativistic medium can be found in Ref. |161 117j . I set this 
issue aside for the present. 

4 Derivation of the vacuum hydrodynamic equations 

The energy-momentum tensor for a perfect fluid has the form [g^ u = diag(l, —1, —1, —1)] 

T nv = + P ) u n u v _ p^w > (4.22) 

where £ and P are respectively the energy density and the pressure in the co-moving 
frame, and is the flow 4-velocity 

^ = ( Uo ,n) = ( 7 , 7 v), (4.23) 

with 7=1/ \/l — v 2 . Energy and momentum conservation equations follow from dJT^ = 
with, respectively, v = and u = i = 1, 2, 3: 

| [{S + P)ul -P} + A J [(S + py «o] = 0, (4.24) 

^ [(£ + P)u^] + A [(£ + PyV + P^] = 0. (4.25) 

Substituting for S and P — from either 1)3.141 13.15|) in the empty- vacuum case, or ()3.171 
13.18)) in the Higgs- vacuum case — yields the fundamental equations of "vacuum hydrody- 
namics." 

Provided the solution is free of shock discontinuities and singularities such as vor- 
tex lines, one may work with equivalent equations that are somewhat neater. (Cf. the 
discussion below Eq. (|2.3j) .) These neater equations arise from u u d^T^ v = ("energy 
conservation in the co-moving frame") and (g au — u a u l ,)d l _ l T lxu = ("Euler equation") 
18 . In the empty vacuum case (in 1 + 1 dimensions) these equations are: 

iip)t + {ipv) x = 0, (4.26) 

((I + gphv) t + ((I + gph) x = 0. (4.27) 

Note that one could re-scale p to p = gp and hence eliminate g from the equations. 

Since one is interested in small disturbances it is natural to assume that v and gp 
are small. I shall refer to this as the nonrelativistic-flow approximation (NFA). In the 
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empty- vacuum case it is equivalent to nonrelativistic reduction. Restoring the factors of c 
and considering 1/c to be a small parameter, e, one obtains the bookkeeping rules for the 
NFA (H] 

v = 0(e), g P = 0{e% ^ = O(e). (4.28) 
Applying these rules to the above equations yields 

Pt + (pv) x = 0, (4.29) 

v t + {\v 2 +g P ) x = 0. (4.30) 

Comparing these equations with (|2.1[) . (|2.3|) . one sees that they correspond to an otherwise 
"normal" medium with a pressure function P(p) = \gp 2 - This is a special case of a 
"polytropic gas" (P oc p n ) with polytropic index n = 2 [U03IIEEI1- Ordinary gases are, to 
a good approximation, polytropic, but with n between 1 and 5/3. However, the n = 2 case 
arises in the treatment of shallow water waves (wavelengths long compared to the depth 
of the water) , where p is proportional to the vertical displacement of the water surface jS] . 

More relevantly, these equations also arise in the description of the free expansion 
of an atomic Bose-Einstein condensate when the trapping potential is turned off |211 
I22j . This fact is natural from the present point of view: the clump of atoms, localized 
with a particular density distribution at t = 0, is just a particular disturbance of the 
empty vacuum. Provided that the initial density distribution is sufficiently slowly varying, 
its spreading will naturally be governed by the vacuum hydrodynamic equations. This 
important example is discussed in more detail in Section 6. 

There is an important conceptual subtlety with the NFA. The original equations are 
Lorentz invariant, while the NFA equations are only Galilean invariant (i.e., only approxi- 
mately Lorentz invariant for small boosts). Thus, using the NFA implies restricting oneself 
to a particular Lorentz frame and its "neighbouring" frames. For a normal medium the 
rest frame of the medium obviously plays this role. However, the vacuum is a Lorentz- 
invariant medium; it has no rest frame. The appropriate frame for the NFA is determined 
solely by the initial conditions. If in some frame the NFA conditions are satisfied at t = 
then they will remain satisfied at all later times. One may trivially take a NFA solution 
and boost it by a large Lorentz boost to obtain an approximate solution to the original 
relativistic equations in which v is everywhere close to 1. Only when the range of v values 
is a significant fraction of unity is it necessary to abandon the NFA and return to the 
relativistic equations, (|4.26l I4~27j) . 
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In the Higgs vacuum case the pressure and the energy density exchange roles. As a 
consequence, time and space in some sense exchange roles. In the 1 + 1 dimensional case 
one obtains 

M x + (rt<rv) t = o, (4.3i) 

((1 + Ca) 1V ) x + ((1 + Cah) t = 0, (4.32) 

which are the same as the empty vacuum case with x «-> t (and p — > a, g — > C). (In 
higher dimensions the "space <-> time" interchange is obviously more complicated. See 
Eqs. (GEHl EHJ) below.) 

Using this symmetry one can define a consistent NFA by the bookkeeping rules: 

v = 0(e), Ca = 0(e 2 ), ^ = 0(e). (4.33) 

ox 

and obtain the equations 

0x + (<n>)t = 0, (4.34) 

v x + (lv 2 + Ca) f = 0. (4.35) 

Although the flow velocity is nonrelativistic (v <C 1), disturbances tend to "propagate" 
superluminally, at 1/v. Hence, the NFA here is not a normal nonrelativistic reduction. 
The resulting equations are "anti-Galilean" invariant. Consider two Lorentz frames, with 
the primed frame moving at velocity w with respect to the unprimed frame. If w is small, 
O(e), then the Lorentz transformations become approximately (j w = 1/vT— w 2 ) 

m=^(m +w m) , *w (436) 

— -7 [ — + w— ) ~ — + w— (4 37) 

dx' w \ dx dt J dx dt 

This is backwards with respect to the usual Galilean transformation, where the x derivative 
would not change and the t derivative would change by aw term. 

This is certainly strange, and takes some getting used to, but one should simply view it 
as an approximation to the full Lorentz transformations, valid in the stated context. One 
is used to dealing with small objects that move slowly, so that their density distributions 
vary rapidly in space, but slowly in time. In the present case one is dealing with large 
objects, slowly varying in space, but relatively rapidly varying in time. This is related to 
the fact that the Higgs vacuum, as a spontaneous Bose-Einstein condensate, has almost 
all its particles in the same quantum state. Small disturbances of this state involve vast 
numbers of particles, spread over long distances, all moving nearly in lockstep, so that 
the disturbance varies only slowly with position while the whole collective has the same, 
relatively rapid time dependence. 
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5 Shock speeds 



Returning to the fundamental, relativistic energy and momentum conservation equations 
one can write down the two equations governing shock discontinuities. In the empty- 
vacuum case one finds 

shock speed = „, M ' +f^ ]] = M^pf^^, (5 . 3 8) 
[i(p(l + gp)^-±gp2]} Ml +gph 2 v]] V ' 

which in the NFA indeed reduces to the usual, nonrelativistic result (j2.11|) with P(p) = 



\9P 2 - 



In the Higgs-vacuum case the corresponding equations are 

allock speed = M^cM = (5 . 39) 

Ml + CVh 2 - all + \Ctr)}} Ha(l+Ca)iH]} ' 

Since ^y 2 v 2 = 7 2 — 1, one may re- write these equations as 

shocks „ - m = BZP+g^fj^l. (5.40) 

Comparing with (|5.38f) one finds that the shock speeds in the two cases are reciprocals of 
each other, in keeping with the x «-> t symmetry. In the NFA the equations reduce to 

shock speed = rr 911 = T &. (5.41) 

6 An exact solution 

In 1 + 1 dimensions, the empty- vacuum hydrodynamic equations in the NFA, (|4.29|l4~3*U|) . 
are the well-known "shallow-water-wave equations," a special case of the polytropic gas. 
Thus, a general, exact 1 + 1-dimensional solution can, in principle, be obtained using the 
method described in Sect. 98 of Ref. [IJ. I shall not pursue that program here, but rather 
I shall discuss a specific 1 + 1 dimensional solution that has a clear physical interpretation 
in terms of the free expansion of an atomic Bose-Einstein condensate — a phenomenon 
has been studied both theoretically [22| and experimentally [23,. From this solution, by 
x t, one can then obtain a solution to the 1+ 1-dimensional Higgs-vacuum hydrodynamic 
equations; a solution that can then be generalized to a spherically symmetric solution in 
d+1 dimensions. 

Ref. |22| has found an exact, analytic solution of the empty-vacuum hydrodynamic 
equations in arbitrary dimensions. Let us consider that solution in 1 + 1 dimensions |24j : 



nn — M 2 l ( h 2 x 2 

\r b—1 X 



< x < b£, (6.42) 
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where M is a normalization constant and b is a function of t given by 



y/b(b-l) + ln(Vb + V&^l) = -jt. 



(6.43) 



It is straightforward to show that b is a solution to the differential equations 



(lb 
dt 




M b-1 



(6.44) 



and 



fb 
dt 2 



2b 2 e 2 



(6.45) 



and hence to verify that Eqs. (j4.29irOU|) are exactly satisfied. Note that the normalization 
constant N should be small and that the scaling of v,gp, d/dt with M is then completely 
consistent with the requirements of the NFA, Eq. (|4.28j) . 

The function b{t) starts from b = 1 at t = and grows linearly for large t. At t = the 
initial flow velocity is zero everywhere and the initial density distribution is an inverted 
parabola between x = — i and x = +£ and zero outside. The solution is self-similar, in 
that the density maintains an inverted parabolic form as its spreads out. See Fig. 3. 

The initial density distribution is in fact the equilibrium density distribution in a 
harmonic trapping potential, in the Thomas-Fermi approximation. Thus, the solution 
with these initial conditions is naturally realized in atom-trap experiments, simply by 
abruptly turning off a pre-existing trapping potential [2~2l 12~3*] . 

In Fig. 3. I have taken v = outside the range — bi < x < +b£. This choice 
is somewhat arbitrary because, of course, v is physically meaningless anywhere where 
p = 0. The fundamental equations are mass and momentum conservation, (|2.1|) and (|2.2|) 
with P = \gp 2 , which are identically satisfied if p = 0, whatever v may be doing. The 
discontinuity in the derivative of p at x = ±M, whether or not one wishes to regard it as a 
"shock," is certainly a place where the hydrodynamic assumption of slow variation breaks 
down. Corrections to hydrodynamics come into play to smooth out these "corners" in the 
density distribution [23]. I discuss these corrections in the next section. 

By making the transformation x <-> t one can obtain an exact solution to the Higgs- 
vacuum hydrodynamic equations. First note that the solution in Fig. 3 can be joined 
smoothly on to its time-reversed solution at negative times. One can then take slices of 
Fig. 3 at various values of x, which becomes the new t. This leads to Fig. 4 which shows 
the time evolution of a a disturbance initially concentrated around x = 0, though with 
long tails out to infinity on each side, and with v = everywhere. The disturbance spreads 
initially somewhat like the empty vacuum case, but then develops a dip at the centre; it 
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Figure 3: Density and flow velocity as functions of x and t for the solution in Eq. (|6.42[) . 
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then splits into two disturbances moving rapidly (superluminally) apart from each other, 
with each one spreading out. Between them is a stretch of perfectly restored vacuum, 
a = 0. 

This solution can be generalized to d + 1 dimensions. The vacuum hydrodynamic 
equations (|4,341 I4.35|) generalize to 

Va + (av) t = 0, (6.46) 

vv t + V.v + a t = 0. (6.47) 
Assuming spherical symmetry, so that v is radial, gives 

a r + (av)t = 0, (6.48) 

v r + {d-l)- + (\v 2 + Ca) t = 0. (6.49) 
r 

Note that the extra term, with respect to the 1+1 dimensional case, occurs in the second, 
not the first, equation. One may now proceed to look for a self-similar solution of the form 



Ca = ^ (V 



br <t< br, (6.50) 



v=Nf 1 -, 

where b and / are functions only of f = — r. Substituting this Ansatz in the first equation 
produces a solution if 

b' 

f = V (6-51) 



db 

d U = 

6(f): 



where b' = Then the second equation reduces to an ordinary differential equation for 

b" +(d- l)£-il=0. (6.52) 

In the 1 + 1 dimensional case this equation, multiplied through by 6', integrates up imme- 
diately to yield (|6.43|) with t — > r and i — > r. In higher dimensions the equation is less 
tractable, analytically. However, it is easy to find a special, power-law solution valid for 




4(d - (6 ' 53> 
With numerical methods one can find solutions that start from b = 1 at r = and 
approach this behaviour at large f. These solutions provide the analog to the 1 + 1- 
dimensional solution discussed earlier, except that the large-r dependence is r 2 / 3 rather 
than linear. The small-f behaviour is 

b = l + i-r 2 + 0(f 4 ). (6.54) 
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a t = 



v t = 




Figure 4: Solution to the 1-dimensional Higgs-vacuum hydrodynamic equations obtained 
by x <-> t from Fig. 3. Only positive x is shown; a is symmetric and v is antisymmetric in 
x. Time is in units of i/N. A qualitatively similar spherically symmetric solution exists 
in higher dimensions. 
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The subleading behaviour at large-f is rather curious. Writing b = bo + (3 and assuming 
-C bo yields a linear equation for (3: 

8' 4 4 B 
P" + (d-l)^ + -(d--)^ = 0. (6.55) 
r 6 3 r z 

Putting r = e z yields an equation with constant coefficients: 

+ (d-2)/3+|(d-|)0 = O, (6.56) 

with solutions (3 = e pz with 

p = - (q - 3d ± a/100 - 84d + 9d 2 ) . (6.57) 

For d = 8 and above the roots are real, but for d = 2, ... 7 the roots are complex. For 
d = 3 one finds 

/?= -^Lcos(4I In r + c 2 ), (6.58) 

where ci,C2 are constants. The presence of the oscillatory factor is, however, all but 
invisible in a plot of b(f). See Fig. 5. Since b is qualitatively similar to the 1+1 dimensional 
case (the main difference being the r 2 / 3 rather than linear growth at large r) the solution 
for a and v as a function of radius at various times is qualitatively similar to Fig. 4. 

3.5 
3 
2.5 

b(f) 2 
1.5 

1 

0.5 

2 4 6 8 10 




Figure 5: The solution to Eq. (16321) with d = 3, for 6(0) = 1, b' = 0. (r = (M/r)r). 
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7 Corrections to empty-vacuum hydrodynamics 



In the atomic Bose-Einstein condensate literature the empty-vacuum hydrodynamic equa- 
tions arise as an approximation to the Gross-Pitaevskii (GP) equation: 



Bib ft 2 

ih-^ = V 2 ip + V(rU + U \ip\ 2 ip. (7.59) 

at 2m 

The constant Uq = Airati 2 /m, where a is the s-wave scattering length is related to the 
constant g introduced earlier by g = Uo/m 2 . The GP equation itself is valid provided 
the medium is dilute; na 3 < 1. In the present context there is no external potential, so 
V(r) = 0, and for simplicity I consider only the 1 + 1 dimensional case. 
Substituting 

l p = JP- e iv W ith v = —<p x , (7.60) 
V m m 

and separating real and imaginary parts, one obtains two equations [251 I26| . One is mass 

conservation: 

Pt + {pv)x = 0, (7.61) 
and the other, after taking d/dx, becomes 

Vt+(\v 2 +gp-A( P ^-\^X) =0, (7.62) 



2 \ p 2 p z 

with A = These two equations reduce to the hydrodynamic ones when the density 

is so slowly varying that one may neglect the "quantum pressure" term proportional to 
A. 

Note that the corrections to empty- vacuum hydrodynamics involve third derivatives, 
rather than the second derivative terms characteristic of viscosity. The corrections are 
non-dissipative and preserve the time-reversal (t — > —t, v — > —v) and parity (x — * —x, 
v — > —v) symmetries of the hydrodynamic equations. 

The two real equations (|7.61j) and (|7.62j) are equivalent to the single complex GP equa- 
tion. It is perhaps helpful to think of this system of equations as the senior member of 
a family of four equation systems, the others corresponding special cases with g or A or 
both set to zero; see Table 1. With g = A = we have the Burgers system, corresponding 
to a pressureless gas of free, classical particles. This system has been studied recently 
by Choquard [27j, who points out the solution p oc v x [2S1- Inclusion of particle interac- 
tions (non-zero g) leads to the empty-vacuum hydrodynamic equations. Returning to free 
particles (g = 0) but including the quantum term (non-zero A) corresponds to the free 
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Burgers 

v t + (^) x = 
Pt + (p v ) x = 


Empty vacuum hydrodynamics 

v t + {±v 2 + gp) x = 
Pt + ((»>)„. = 


Schrodinger 

Pt + (pv) x = 


Gross-Pitaevskii 

Vt + (V +gp-A(£f-\ji)) x = Q 

Pt + (pv) x = 



Table 1: Family of four equation systems related to empty- vacuum hydrodynamics. A = 
h 2 /{4m 2 ). 

Schrodinger equation, which is manifestly linear in its usual complex form. With both g 
and A non-zero one has the nonlinear GP equation. Appendix A discusses the infinite 
tower of conservation laws associated with these equation systems. 

8 Corrections to Higgs vacuum hydrodynamics 

What are the appropriate higher-derivative corrections to the hydrodynamic equations in 
the Higgs-vacuum case? To answer this question, in the absence of a microscopic theory, 
necessarily involves some guesswork. Nevertheless some simple physical considerations, 
together with a mathematical/aesthetic criterion appear to point to a unique choice. 

The physical considerations from the analogy with the empty vacuum case are (i) 
one expects third-derivative, non-dissipative, time-reversal invariant terms; (ii) in the 
NFA these will be time derivatives because space derivatives count as O(e), (hi) the new 
constant multiplying the new terms should count as 0(e 2 ). However, there is a further 
physical consideration that suggests that the new terms should break the x <-> t relation 
to the empty vacuum case: unlike the empty-vacuum case, v is meaningful even when 
a = 0. That is because in the Higgs case there is a pre-existing density of particles: a 
disturbance that causes these particles to move non-uniformly will cost energy, even if it 
nowhere changes the density at some initial instant. 

Guided by these considerations and a mathematical/aesthetic criterion that the new 
equations, in 1+1 dimensions, should be an integrable system, I looked for equation systems 
that would generate a tower of conservation laws. In this way I was led to the following 
equations: 

v x + {±v 2 + Ca) t = 0, (8.63) 
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a x + (av - Tlv u ) t = 0. (8.64) 

The second of these equations represents momentum conservation. As discussed in Ap- 
pendix A, these two equations imply other conservation laws. I conjecture that there are 
an infinite number. The next of these, in its simplest form, is 

(av) x + (av 2 + \Ca 2 + ll(\v 2 - vv tt )) t = 0. (8.65) 

Alternative forms can be obtained by adding a total x derivative to the density and 
subtracting the corresponding t derivative from the flux. In this way one can obtain an 
energy-conservation equation whose energy flux coincides with the momentum density: 

(av - Kv u ) x + (av 2 + \Ca 2 - Tl(\v 2 + 2vv tt + Ca tt )) t = 0. (8.66) 

It appears that the constant 1Z must be positive, so that the new vut term in Q8.64JI 
has the effect of smoothing out the discontinuity in v and the corresponding "corner" in 
a in the solution shown in Fig. 4 |29j . 

Now, in analogy with Table 1, one has a family of four equation systems; see Table 
2. The pattern is somewhat different from Table 1, where the mass-conservation equation 
Pt + (pv) x = held in each case. Its x <-> t analog in Table 2 (which is not mass conserva- 
tion, but momentum conservation) becomes modified in the higher-derivative cases in the 
bottom half of the table. 



RIF 

+ (lA = 
a x + (av) t = 


Higgs-vacuum hydrodynamics 

v x + (±v 2 + Ca) f = 
a x + (av) t = 


Schrodinger analog 

v x +(\v 2 ) t = Q 

Ox + (o-v - Kvtt) t = 


GP analog 

v x + (iv 2 + Ca) t = 

(Tx + ((TV - Hv U )t = 



Table 2: Family of four equation-systems related to Higgs-vacuum hydrodynamics. 

The simplest case (top left in Table 2) is denoted "RIF" for "relativistic incompressible 
fluid" since it corresponds to zero compressibility, C = 0. Its Euler equation v x + vvt = 
is d^u^ 1 = 0, the relativistic version of the usual incompressible-fluid condition V.v = 
[I]. The higher derivative cases, for lack of better names, are denoted as "Schrodinger 
analog" and "GP analog," even though the nature of the analogy is not clear. 
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I suspect that in the "Schrodinger analog" case there is some function VP((T, cp), with 
(ft = v, that satisfies some linear equation, but I have been unable to find this function 
and equation. I have found, though, that 

' = -?®, ^ 

yields a solution to the a equation of the "Schrodinger analog" system. More generally 
one can add to a a term proportional to v%, since this is a solution in the RIF case. (Cf. 

nam.) 

Introducing a velocity potential, ip, with ipt = v, in the "GP analog" system allows 
one to solve the first equation for <r: 

Ca = -(<p x + ^ 2 t)- (8-68) 
One may then use this to eliminate a in the second equation to obtain 

3 2 



ip xx + 2(p t (p tx + tp tt i-p t + tp x J + KC<putt = 0. (8.69) 
Note that only the product of 1Z and C appears in this equation. 



9 'Bright' and 'dark' solitons 

For the GP equation it is well known that in the negative g case, corresponding to attractive 
particle interactions, there are "bright soliton" solutions where the density has a sech 2 
form |25| . (Such solitons have been observed experimentally |3Uj.^ Remarkably, there are 
similar solutions to the "GP analog" equations, (|8.63l 18.64(1 if 1ZC is negative. Although 
this "wrong sign" case probably has no physical significance these solutions are of 
some interest because (i) they give support to the conjecture of integrability, and (ii) they 
hint at a symmetry, more subtle than x <-> t, between the GP and "GP analog" equations. 
To begin, one postulates a travelling-wave solution 

v = v(x — at), a = a(x — at), (9.70) 

where the constant a represents the soliton's velocity. Substituting in Eq. (|8.63|) and 
integrating with respect to the variable s = x — at gives 

(1 - f 2 ) 

v - a(\v 2 + Co) = const. = v ; . (9.71) 
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The constant of integration has been written in this form for later convenience. Note that 
when (7^0 the velocity v tends to one or other of the constant values v± = (1 ± £)/a. 
Substituting Eq. (j9~7U|) into Eq. (|8~M|) one obtains 

cr{l - av) = -Ka 3 v", (9.72) 

where the constant of integration must be zero to have v" = when a = 0. Substituting 
for a from the previous equation leads to 

(£ 2 - (1 - avf){l - av) = -2UCa b v". (9.73) 

Multiplying through by av' and integrating yields 

4KCa 6 v' 2 = av(2 - av)(2(l - £ 2 ) - 2av + aV) + const.. (9.74) 

The constant of integration here must be —(1 — £ 2 ) 2 so that v' —* when v — > v± = 
(1 db Cj/a. The result then simplifies to 



v 



'2 _ __}__ ( (~\ „„.\1 c2\ 2 

AKCd 



{(1-avf-ef- (9-75) 



Since the left-hand side is positive, a solution of this type is possible only if 1ZC is negative. 
It is convenient to define 



k = V-iKC. (9.76) 
Taking the square root of Eq. 1)9. 75j) and integrating yields the solution 



1 ± £ tanh ( — ^(x — at) 
1 Ka z 



(9.77) 



From Eq. (|9.71|) one obtains the corresponding form of a: 



; r sech 2 (-^(x -at)) . (9.78) 



2Ca 2 \Ka 

The above solution is characterized by two parameters £ and a and by the sign choice, 
db, in Eq. (|9.77[) . Each requires a little discussion, (i) Without loss of generality £ may 
be taken as positive. For £ = the solution becomes trivial (a = and v = constant, 
representing the vacuum itself). For a single soliton one could set £ = 1 by choosing 
a frame of reference in which V- = 0. However, for the multi-soliton case, discussed 
below, the extra generality of the £ parameter is important, (ii) The a parameter can 
have either sign but, for the NFA to be valid, the flow velocity must be small everywhere 
and hence 1/ | a \ <C 1. Hence, the soliton moves at a hugely superluminal speed | a 1. 
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The total momentum and energy carried by the soliton are, however, small and inversely 
proportional to its speed and speed-squared, respectively: 

^soliton = M f~) , ^soliton = M (~) ' ( 9J9 ) 

with 

f°° K 

M= dxu = i—. (9.80) 

J-oo ^ 

These results follow from evaluating the momentum and energy densities, which turn out 
to have the same sech 2 form as a; that is, using Eqs. ()9.77l l9~?8*)) . one finds that 

av — IZvu = cr/a, (9.81) 

av 2 + \Ca 2 - Il(\v1 + 2vv tt + Ca tt ) = a/a 2 . (9.82) 
(iii) Another important property of the solution is that 

° = T^v t , (9.83) 

with the upper /lower sign determined by the upper/lower sign in ()9.77|) . I shall refer to 
the lower-sign choice as "solitons" and the upper-sign choice as "anti-solitons." Solitons 
have positive vt while anti-solitons have negative vt- 

Multi-soliton solutions can be constructed by using the ansatz a = -^vt- Both of the 
"GP analog" equations are then satisfied if 

v x + vv t = --v u , (9.84) 

which is the viscous Burgers' equation (|2.12l) with x <-> t. This equation can be solved by 
the Hopf-Cole transformation |l()j . One defines 

i) = e^, with (p t = v, (9.85) 

and requires ip to satisfy the linear equation 

V'x = ~\^tt (9.86) 

(the diffusion equation with x <-> t). The resulting v will then satisfy (|9.84j) . The single- 
soliton solution, ()9.771 19.78|) with the lower-sign choice, satisfies the a = -^jVt ansatz and 
by integrating v one can find its corresponding ip. The appropriate constant of integration 
is a function of x chosen so that ip = e v ^ K satisfies Eq. (|9.86|) . Hence, one obtains the 
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t/j function corresponding to a single bright soliton. By superposition, one can then write 
down a general class of solutions: 

ip = V]i^ n exp ( -^—] exp f-(x - gn) ^" 1 " n"^ cosh - x n - a n t)] , (9.87) 

where K n , a n , £ n , and x n are arbitrary constants. By differentiation one can then obtain 
v and a as 



"="f- '=M?A' (a88> 

Multi-anti-soliton solutions can be constructed in the same way by reversing the sign 
of k in the above. However, it is not clear if one can analytically construct solutions 
involving both solitons and anti-solitons. 

As noted earlier, these 'bright soliton' solutions are probably unphysical since they 
require TZC < In the physical case, 1ZC > 0, there are analogs of the 'dark soliton' 
solutions to the positive-^ GP equation. One may find these by the same method as 
before, but in Eq. Q9.74JI one now needs a different constant of integration. Choosing it 
to be 2£ 2 — 1 yields a solution where v tends asymptotically to 1/a at ±oo. The resulting 

solution is 

' / ' (9.89) 



1 

v = — 
a 



a = a 



1 ± \J Aa 2 Cao sech — a ^^j 

1 - 2sech2 (^^" at) ) 



(9.90) 



where o"o = S, 2 / (2a 2 C) is an arbitrary constant. In this solution both v and a are constant 
at ±oo but have a localized "dip" or "bump" that moves with velocity a. The dip/bump 
in v can be arbitrarily small, but a has a dip that always goes negative, with a minimum 
value —<to. (See footnote ^2] for remarks on negative a.) 

Since this 'dark soliton' solution does not have a — > at ±oo it does not represent a 
perturbation of the vacuum. However, one can regard it as describing a localized distur- 
bance in a larger, and much more slowly varying, perturbation of the vacuum. Thus, a 
solution like that shown in Fig. 4, for example, could be "decorated" with one or more 
'dark soliton' dips. Calculating the momentum and energy densities for the dark soliton 
one finds, after subtracting the background values associated with a = a® and v = 1/a, 
that they are (a — o"o)/ a and (a — o"o)/a 2 , respectively. 

Note that the above solution can always be boosted to another frame with an "anti- 
Galilean" transformation to obtain another solution. In particular one can boost to a 
frame where v — > asymptotically. In that frame v and a become functions of t only. 
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10 Discussion 



There are two main points that I wish to emphasize: (i) Hydrodynamics in the empty- 
vacuum case makes perfect sense and describes an experimentally observed phenomenon, 
the free expansion of an atomic Bose-Einstein condensate when the atom-trap potential is 
turned off. (ii) Hydrodynamics in the Higgs-vacuum case gives very strange and exciting 
behaviour as a consequence of the fact that the speed of sound in the Higgs vacuum is 
formally infinite. The Higgs vacuum is a medium that is both ultrarelativistic (pressure 
^> energy density) and ultra-quantum, being a Bose-Einstein condensate with almost all 
its particles in the same quantum state. Not surprisingly, perhaps, its properties are very 
different from those of familiar media. 

The vastly superluminal speed of hydrodynamic disturbances in the Higgs case will 
obviously require careful study. It may be important to point out that this phenomenon 
can occur without any material particle moving faster than the speed of light; see Fig. 
6. In fact, a similar issue arises even for a non-relativistic Bose-Einstein condensate |3 1 j : 
phonons with low momentum, k — > 0, travel at a finite speed of sound, yet they are 
made up of atoms that move with very small velocity k/m. In the non-relativistic case 
a spacetime picture like Fig. 6 predicts that the number of excited atoms present at any 
one time equals the ratio of the speed of sound to the speed of the individual atoms, a 
result that indeed follows from the usual Bogoliubov theory. See Ref. |31j for a detailed 
discussion of this point. 

There is an apparent dichotomy between hydrodynamics, which implies "soft" modes 
that can be excited with arbitrarily little energy, and quantum field theory which (unless 
it explicitly contains massless particles) predicts an energy gap. In the empty-vacuum 
case the atoms carry a conserved quantum number, so there is a superselection rule fixing 
the total number of atoms. The theory splits into separate, non-communicating "sectors," 
each with a different total number of atoms. The vacuum, strictly speaking, is in the 
zero-atom sector and its lowest excitation, an atom/anti-atom pair, requires an energy 
2Mc 2 , where M is the mass of the atom. However, if one is given a system with N atoms 
then the relevant "vacuum" state — the lowest energy state with N atoms — has those 
atoms existing but dispersed in all directions to infinity (assuming a repulsive interaction 
g > 0). This "vacuum" can be perturbed by bringing some or all of these atoms to finite 
distances from each other, and this may be done with arbitrarily little energy. In this 
sense, there are soft, hydrodynamic modes. 

The situation for the Higgs vacuum is different. No conserved quantum number is 
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X 



Figure 6: Possible spacetime picture of a superluminal phonon-like excitation. There are 
multiple pair-creation and pair-annihilation events organized so that the collective motion 
is fast L/r » c while the individual particles move slowly (L/(2T) < c). 



carried by the particles that spontaneously condense (the quanta of the scalar field that 
acquires a vacuum expectation value). According to conventional theory, the only excita- 
tions are massive (quasi)particles, and hence the system has a gap. However, the equation 
determining the inverse propagator at zero momentum actually has two solutions; M?, 
where is the Higgs mass, and zero |32| . This fact suggests that there may be some 
deeper, underlying theory that possesses gap-less, hydrodynamic modes. In the limit 
where this theory reduces to the conventional field theory, the realm of these hydrody- 
namic modes would shrink to p^ 1 = 0. Studying the hydrodynamic equations is then, I 
believe, a window into this unknown, deeper theory. 

An example of what I have in mind is provided by the sine- Gordon equation |33| 

d 2 e d 2 e r 
w -— 2+sm e = o. (10.91) 

For small 6 this equation reduces to the Klein-Gordon equation and hence the dispersion 
relation is u(k) = \f \ + k 2 . One might therefore think that the spectrum of excitations 
necessarily has a gap. However, the non-linear equation possesses "breather" solutions 

9{x,t) = 4arctan ( ^t/VT+^)_ ) , (10 . 92) 

V cosh(e(x — xq)/V 1 + e 2 ) / 



whose fundamental frequency, ujb = + < 1 lies in the gap. (The sine-Gordon 

"breather" is unique in being stable. Similar breathers exist in the (f> case but they decay 
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due to e -1 / e couplings to linear, plane-wave modes [SI]-) The moral is that even if the 
linearized theory has only massive excitations, one cannot discount the possibility of soft, 
non-linear modes. 

In this paper only the classical hydrodynamical equations have been considered. A 
task for the future is to re-cast the equations into Hamiltonian form so that one can apply 
canonical quantization. For a normal medium in the linear, acoustic regime this procedure 
would take one from classical sound waves to phonon quanta. For the vacuum case the 
quantization will be inherently non-linear. 

Acknowledgements I thank Maurizio Consoli for discussions. This work was sup- 
ported in part by the U.S. Department of Energy under Grant No. DE-FG05-97ER41031. 
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Appendix A: Towers of conservation laws 



For each of the 1 + 1 dimensional equation systems in Table 1 there is an infinite tower of 
conservation laws 

D t (n) + f]"' = 0, (10.93) 

where D is a density and F is the associated flux. Each will be normalized such that 
its pv n ~ l term has unit coefficient. For n = 1,2,3 the conserved quantities (J dxD^ 1 ) are 
mass, momentum, and twice the nonrelativistic energy. Of course, one has the freedom 
to add a total x-derivative to a density; i.e., D' = D + T x and F' = F — T t satisfy 
D' t + F' x = 0. 

The Burgers case is trivial, since it corresponds to non-interacting classical particles. 
One has simply 

£>(») = pv n -\ = pv n . (10.94) 

That is, not only are mass, momentum, and (nonrelativistic) energy (integrals of p, pv, 
and \pv 2 ) conserved, but — since v is constant for a free particle — so is the integral 
of pf(v) for any function of v. (Thus, the set of conservation laws above is really only 
a subset of a continuous infinity of conservation laws. I focus on this subset because it 
generalizes to the other equation systems.) 

The empty-vacuum hydrodynamic equations include the effects of particle interactions 
(repulsive for g > 0). The D's and F's are now polynomials in p and v: 

£>(») = v n+1 d n (gp/v 2 )/g, = v n+2 f n (gp/v 2 )/g (10.95) 

where the functions d n (s) and f n (s) (with s = gp/v 2 ) satisfy 

(n + l)d n = f n - (1 - 2s)d' n , (10.96) 

(n + 2)/ n = (1 + 2s)f' n - (1 - s)d' n , (10.97) 
whose polynomial solutions are 

[("+i)/2] / _ n , ! 

"»W= E ( n -a + i)!i!(i-i)! *'' ( la98 > 

1=1 

[(n+2)/2] _ 

/.(.)= e '"+ 1 -') ( „'! a + 2) , f 10 -") 

Two noteworthy relations are 

£,(n) = (10.100) 

n c/f 
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F (n) = l(£,(n+l) + vD {n)y (10.101) 

In the Schrodinger case one goes back to free particles, but now they are quantum 
mechanical. The D : s and F's now involve higher derivatives. They can be expressed 
compactly in terms of the wavefunction ip. Defining if)\m as the iVth derivative ip xx ... x one 
has 

D {2N+l) =r [N] ip[N], (10.102) 

F (2N + 1) = D[ 2N + 2) = Zl ^ [N+1] - r [N+1] M) > (10-103) 

= ~ (2r [N] M - r [N -^[N+i] - r [N+ i^[N-i]) . ao.m 

By substituting if) = p/m e lLp with v = ^<f x one can obtain the -D's and F's in terms 
of p, v, and their derivatives, though the expressions quickly become cumbersome. (See 
Table 3 with g = for the first few.) There is a straightforward generalization of these 
results to d+1 dimensions, with the D^'s being alternately scalar and vector densities and 
the fW's correspondingly vector and tensor fluxes. I have not seen these conservation 
laws for the free Schrodinger equation in any quantum-mechanics textbook, although they 
surely must be known. 

In the Gross-Pitaevskii (GP) case one has both quantum effects and interactions. It 
is a highly non-trivial property of this equation, related to its integrability, that it admits 
an infinite number of conservation laws in the 1 + 1 dimensional case |35| . The first four 
of these are given explicitly in Table 3. 



n 




Fin) 


1 


P 


pv 


2 


pv 


pv 2 + \gp 2 - A [p xx - f) 


3 


pv 2 + gp 2 + A& 


pv 3 + 2gp 2 v + A (&£!L + 2 Px v x - 2p xx v) 


4 


pv 3 + 3gp 2 v 


pv 4 + \gp 2 v 2 + g 2 p 3 




+ A ( P p + 2 PxVx ZPxxV^J 


+ A (-5p xx v 2 + 6^-v 2 + 4p x vv x + 2pv 2 - 2pvv xx \ 






<Ag{if%- 3ppxx) 






i A 2 f-Zpxxxpx i 2plx i Pxxpl p%\ 

+ A { p 1 p 1 p' ^) 



Table 3: The first four conservation laws, Dt + F x = 0, for the GP equation system. 
A = h 2 /{4m 2 ). is the mass density, is the momentum density, and is twice 
the nonrelativistic energy density. 
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I turn now to the equation systems of Table 2. The RIF and Higgs- vacuum hydrody- 
namics cases are related by x <-> t to the Burgers and empty- vacuum cases, respectively 
and so inherit the same set of conservation laws as 

+ F t (n) = 0. (10.105) 

(Note that I choose not to switch the names D and F, so one needs to remember that now 
D is a flux and F is a density) 

For the "Schrodinger analog" case (JZ / 0, but C = 0) it is easily verified that there 
are an infinite number of conservation laws that need only involve a, v, vt, vu- 

F n = av n + \U ((n - \fv n - 2 v\ - 2v n - l v tt ) , (10.106) 

D n = ov n - 1 + \U(n - l)(n - 2)v n - 3 v? (10.107) 

(One may of course modify these by adding a total t-derivative to D and subtracting the 
corresponding x-derivative from F. I have chosen to do this for the n = 2 case in Table 4 
in order for D^ 2 \ interpreted as the energy flux, to equal the momentum density F^\) 
Table 4 shows the first four conservation laws for the "GP analog" system. 



n 






1 


a 


av — TZvu 


2 


av - IZvu 


av 2 + \Ca 2 - K{\v 2 + 2vv tt + Ca tt ) 


3 


av 2 + Co 2 + Tlv 2 


av 3 + 2Ca 2 v + TZ(2vv 2 - v 2 v tt + 2C(a t v t - av tt )) 


4 


av 3 + 3Ca 2 v 


av 4 + \Ca 2 v 2 + C 2 a 3 + Tl{\v 2 v 2 - v 3 v tt ) 




+ K(3vv 2 - ACav tt ) 


+ lZC(6atWt — 3av 2 — lOavvtt) 
+ 2K 2 Cv 2 tt + nC 2 (2a 2 t - 4aa tt ) 



Table 4: The first four conservation laws, D x + F t = 0, for the "GP analog" system 
(Higgs- vacuum hydrodynamics plus higher-derivative corrections). is the momentum 
density and F^ is the energy density. 
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